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A model is formulated to predict numerically the pressure 
distribution induced by a round, turbulent, unheated, subsonic 
jet exhausting normally through a flat plate into a subsonic 
crossflow. The complete model assumes that the predominant 
features of the flow are jet entrainment and a pair of contra- 
rotating vortices which form downstream of the jet. Experimentally 
determined vortex properties and a reasonable assumption concerning 
jet entrainment are used. Potential flow considerations are 
used ^cept in the wake region, where a single method for 
approximating the pressure distribution is suggested. The 
calculated pressure distribution, lift, and pitching moments 
on the flat plate are presented for a jet to crossflow velocity 
ratio of 8 and are compared with experimental results. 
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INTRODUCTION 


During the transitipn from hovering to horizontal flight, 

a VTOL aircraft may experience loss of lift and undesirable 

; 

pitching moments due to the effects of one or more jots exhausting 
normally into the crossflow created by the craft’s forward motion. 
Figure 1 is a sketch which gives a qualitative description of the 
flow field resulting from a single- jet VTOL aircraft transitioning 
from hovering to horizontal flight. The interaction of the jet 
and crossflow induces a pressure distribution on the lower surfaces 
of the craft, resulting in loss of lift and, in most cases, a 
nose-up pitching moment. Although most VTOL craft employ multiple 
jets which are of a higher ten^erature than the surrovinding air, 
it is. generally desirable to sinplify the problem to a case 
which lends itself to experimental and theoretical analysis 
while still retaining the basic characteristics of the jet-and 
crossflow interaction process. Therefore, it is convenient 
to restrict the problem to that of a single, unheated, subsonic 
jet exhausting at large angles through a flat plate into a subsonic 
crossflow. This simplified case has applications in problems 
other than the transitioning of VTOL craft, such as the problem 
of cooling gases in turbine combustors and the discharge of 
effluent or cooling water into a waterway. 

Early investigations into the problem of a jet in a cross- 
flow were concerned with studies of the jet deflection resulting 
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from the interaction of the Jet and crossflow. These studies 
resulted in empirical relationships for the Jet centerline which 
is defined as the locus of points of maximum velocity in the 
symmetry plane (ref. 1). Recent investigations into the velocity 
field associated with a Jet in a crossflow have shown the Jet wake 
region to be dominated by a pair of contrarotating vortices which 
form downstream of the Jet and to either side of the symmetry 
plane. Fearn and Weston (ref. 2) investigated the case of an 
imheated Jet exhausting normally into a crossflow and formulated 
empirical relationships for the characteristics of the vortex 
pair , 

The pressure distribution on the flat plate lias received 
considerable attention by several investigators. Fearn and 
Weston (ref . 3) , Thompson (ref . 4) , and Bradbury and Wood (ref . 5) 
investigated the flat plate pressure distribution for several jet 
to crossflow velocity ratios. Wooler (ref, 6) and Thompson each 
attempted to predict numerically the fiat plate pressure distri- 
bution by modelling the various features of the flow with potential 
flow singularities. The development of potential flow models was' 
limited by the lack of data concerning the characteristics of the 
vortex pair. Also, the use of potential flow modelling alone resulted 
in erroneous pressure calculations in the wake region downstream 
of the Jet. 

This paper will present a method for calculating the flat 
plate pressure distribution associated with a Jet in a crossflow. 

As in the models presented by Thompson and by Wooler, potential 
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flow singularities are incorporated. However, in tMs study, the 
model incorporates the physical characteristics of the vortex 
pair, as described by Fearn and Weston (ref. 2) and extended 
by Sellers (ref. 7). As a result, the model is more consistent 
with the flow field actually observed than were the previous 
modelling attempts. In addition, a means of approximating 
the pressure distribution in the wake region is suggested. 

Therefore, a fairly complete descri^ition of the flat plate pressure 
distribution, based on the physical characteristics of a jet in a 
crossflow, may be accomplished. The model is applicable for all 
velocity ratios, as long as a graphical or analytic description 
of the vortex properties is available. For this study, it is desired 
to demonstrate the applicability of this model when used in con- 
junction with an analytic description of the vortex pair. Since 
an adequate analytic description of the vortex properties is 
available only for a jet to crossflow velocity ratio of 8 (ref. 7), 
the model presented in this thesis is restricted to this case. 



LITERATURE REVIEW 


In order to fully appreciate the problems encountered in 
the analysis of the jet in a crossflow, it is necessary to review 
briefly some past attempts at modelling the flow and predicting 
the flat plate pressure distribution. Generally, attempts to 
predict the pressure distribution on the flat plate have been 
limited in success, due to the lack of data concerning entrainment 
by the jet and the nature of the vortex pair. 

Numerical predictions of the flat plate pressure distribution 
generally rely on the use of potential flow singularities to 
model the jet and vortex pair. The pressure field is related 
to the numerically calculated potential or velocity field through 
JBernoulli's equation. Generally, no account can be made for the 
wake region throi^h the use of potential flow modelling alone, 
since energy losses throi^h viscous effects and separation result 
in the inapplicability of a potential flow solution, 

Thompson (ref. 4) attempted to predict the flat plate 
pressure distribution by modelling the jet and vortices with a 
distribution of sinks and doublets placed along the jet centerline. 
Since the vortex spacing is generally small compared to the distance 
to a typical field point, Thompson assumed that the vortex pair could 
be modelled as a distribution of finite horseshoe vortices distributed 
along the jet centerline. Since such a system of vortices is 
equivalent to a distribution of axially nca*mal source-sink doublets 
(ref. 4 and ref. 8, pg. 96), Thompson placed a distribution of 
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doublets along the jet centerline. To account for mass entrainment 
by the jet, a distribution of line sinks was also placed along the 
jet centerline. Although Thompson recognized that the vortex 
trajectories lie to either side of the symmetry plane and that their 
projection onto the symmetry plane does not coincide with the jet 
centerline, the doublets were placed along the jet centerline for 
analytical convenience. The doublet strengths were obtained through 
experimental measurements of the velocity field, while the strength 
of the sink distribution was assumed constant and its value 
determined by matching experimental and predicted upstream pressure 
contours. Thompson analyzed several cases, for velocity ratios of 2, 
4, and 8, both with and without the sink distribution. The pressure 
distribution was accounted for fairly well by the model incorporating 
both the sink and- doublet distribution. The model incorporating only 
the doublet distribution could not account for the pressure 
distribution upstream of the jet orifice. As expected, both models 
failed in the wake region. 

Wooler (ref. 6) attempted a numerical analysis of the flat plate 
pressure distribution using two different potential flow models. 

The first, which Wooler called the vorticity model, idealized the 
jet as a system of horseshoe vortices placed along the jet center- 
line. The strength of each vortex was determined from the jet 
momentum flux and the radius of curvature of the jet centerline. A 
contour plot of the pressure distribution on the flat plate shows 
good agreement with the experimental results of Bradbury and Wood 
(ref. 5) except for the wake region and the area upstream of the 
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jet orifice. 

Wooler's second method, called the sink-doublet model, utilized 
a system of sinks and doublets distributed along the jet centerline. 
The doublets institute a blockage of the free stream by the jet, 
essentially modelling the jet as a cylinder of elliptic cross 
section. The sink distribution along the jet centerline models 
entrainment of the surrounding air by the jet. The sink strengths 
were determined by solving equations for the jet development such 
that they provide a best fit with experimental test data. The 
doublet strengths were obtained from the complex velocity potential 
for two-dimensional flow past an ellipse. Wooler did not attempt 
to calculate the pressure contovirs on a flat plate using this model, 
but instead calculated wing loadings caused by a jet issuing from 
the center of a wing. The numerical results for particular 
spanwise and chordwise data points correlate well with Wooler's 
experimental data. However, the choice of data points .did not 
include the wake region or the 'region upstream of the jet, ' 

In the models presented by Wooler and by Thompson, singularities 
modelling the vortex pair were placed along the jet centerline 
rather than along the true vortex trajectories. In addition, 
vortex strengths were, in the case of Wooler's model, obtained 
from assumptions concerning growth of the jet plume rather than 
from an experimental description of the vortex pair, 

A recent experimental study by Fearn and Weston (ref. 2) has 
provided a quantita.tive description of the vorticity associated 
with a jet in a crossflow. This information was found to be 
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useful in the formulation of a more realistic model for calculating 
the pressure distribution on the flat plate. The authors presented 
two models for the vortex pair. The first, known as the filament 
model, idealized the vortex pair as two filament vortices. It 
was assumed that the characteristics of the vortex pair change 
slowly along the vortex curve. Hence, the velocity field at a 
cross section perpendicular to the vortrac curve was assumed to be 
equivalent to that produced by two infinite filament vortices 
perpendicular to the cross sectional plane. From experimentally 
determined velocities in the symmetry plane, the strength and 
spacing of the filament vortices at the cross section were obtained. 
The second model, known as the diffuse model, relaxed the assumption 
of concentrated vorticity and assumed the vorticity distribution of 
each vortex to be Gaussian in nature. Figure 2 shows the coordinate 
system used in the description of the diffuse model. The.vorticity 
at any point was assumed to be 






(i) 


where o> is the maximum vorticity in the cross section, and a 
o 

is the diffusivity of the vortices. The integrated strength 
of a single diffuse vortex was defined as 

Atf CO 


n, = 



-ar 

Uo& rJrJe 


(2) 


'o /c 

The net effective strength P of each diffuse vortex was defined 
.as the net flvix of vorticity across the half plane of the cross 
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section 


1% 



P = j i cj(r,0) rdrJe 


<3) 


- 77/2 > 

The effective spacing h. was defined to be the center of vorticity 

y' j 

I 


k-f^j j Yy rdrdd ( 4 ) 

Equations (3) and (4) were evaluated using the assumed distribution 
of vorticity to give 


r = C erf(pK) 


(5) 


and 


h - \\o/ etfi^ho) 


( 6 ) 


where 

erf(pkc) e Jt 

O 

is the error function.- The authors found that the parameters of 
both vortex models were functions primarily of the jet to crossflow 

I r 

velocity ratio, R, and of arc distance along the vortex curve, S/D. 
The results of .this investigation were presented graphically. It 
was observed that the vortex strength determined from the diffuse 
vortex model was not a function of S/D and could be written as 
a linear function of R, The expression for the integrated 
strength was written 




m 
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where A is a constant equal to 0.72 and ^ is the integrated 
strength, non- di mens ionali zed by 2DU^ for convenience. The 
quantity 2DU^ was found by Chang-Lu (ref. 9) to bo equal to the 
roll up of the vorticity around a two dimensional cylinder. 

Equations (5), (6), and (7) represent a desorption of the strength 
and spacing of the vortex pair, but contain two parameters, p, and h^, 
which must be specified for the description to be complete. 

In this study, it is desired to demonstrate the use of the model 
in conjunction with analytic descriptions of the vortex pair , even 
though the model may be used with either graphical of analytic 
descriptions of the vortices. At this time, analytic relationships 
for ^ and h^ have been developed only for a velocity ratio of 8, 

Sellers (ref. 7), in an extension of the work by Pearn and Weston, 
formulated relationships for ^ and h^ as functions of s/D and R=8. 

The equations may be written 

and 

^ = c(l~ e’) C9) 

where B=2.11 and C=2.04 . With these analytic descriptions 
for p and h^, the effective strength P and spacing h of the 
vortex pair for a velocity ratio of 8 may be calculated, as a 
funttion of s/D. 

Equations (8) and (9) are based. on the assumption that the 
vortices begin at the center of the jet orifice. Since most of 
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the velocity field data was taken at S/D>5, there is considerable 
doubt as to the nature of the vortex trajectories close to the 
jet orifice. Sellers also presented an equation for h^ which 
would result in the vortices intersecting the flat plate at the 
edges of the jet orifice. The equation, again valid only for 
R=8, may be written 

^ = c'(/ - e'^) + •■5' <“> 

where C’ =1.389 . 

Fearn and Weston (ref. 2) also formulated relationships 
for the jet centerline and vortex curve. An equation of the form 

Z/D^ ( 11 ) 

describes the jet centerline and the vertex curve adequately. For 

the jet centerline, a =0.9772, b =0.9113, and c =0.3346 . For the 
^ c c c 

vortex curve, a =0.3473, b =1.127, and c =0.4291 . • 

The equations formulated by Fearn, Weston, and Sellers 
give a fairly complete analytic description of the vortex pair 
for R=8, providing a basis for the formulation of a realistic 
model for calculating the flat plate pressure distribution induced 
by a jet in a crossflow. 



FORMUIATION OP THE PROBIEM 


TH© relationships for the characteristics of the vortex 
pair introduced hy Fearn and Weston (ref, 2) and extended by 
Sellers (ref. 7) make possible the formulation of a model for 
the vortex pair which may be used in a numerical analysis conducted 
through the use of a digital coaiputer. The relationships are 
repeated for convenience. 


' r = r, erHph.] 


r,= 2DU^AR- 2ADUj 'w 

Z!D = af?‘(x/Dr ■ <“> 

~ (s/p)''‘ . 

Equations <5) , (6) , <7) , and (11) are general in that they are 
valid for all s/D and B within the range of the experiment, that 
is,- for 5<s/P<45 and 2<B<10, Equations (8) and (9) are, as pointed 
out earlier, valid only for B=8. 

The vortices are modelled in the cosputer program by a 
distribution of, finite filament vertices placed along the vortex 
trajectories described by equations (6) and (11) . The strengths of 
the finite filaments are given by equation (5). Although it is well' 
'established that the vortices are diffuse in the actual case of a 
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jet in a crossflow, it is assumed that the velocity induced at 
a field point by a filament vortex will he indistinguishable 
from the velocity induced by a diffuse vortex, as long as the 
distance from the filament is large compared to the radius of 
the diffuse core. ^ Since the diffuse vortices are observed to 
vary in strength along the vortex curve as a result of diffusion of 
vorticity across the symmetry plane, the strengths of the filament 
vortices are varied as a function of s/D in the computer model. 
Although this is inconsistent with Helmholtz's Laws, it should 
be emphasized that the filament vortices are used for analytical 
convenience and are modelling a diffuse vortex pair, for which 
there are no restrictions regarding the variation in strength 

along the vortex trajectories. 

The distribution of filament vortices is composed of finite 
straight line segments, the lengths of which are determined by 
the local radius of curvature of the vortex curve. Enough vortex 
segments are generated such that the addition of more segments 
results in negligible velocity changes at the plane of the flat 
plate. ' 

A distribution of line sinks is placed along the experimentally 
determined jet' centerline, defined by equation (11). Although 
little data is available concerning entrainment by a jet in a 
crossflow, experimental and theoretical analyses are available for 
a free jet which may provide a basis for estimating the amount of 
entrainment by a deflected jet and thereby establishing the 
strengths of the line sinks. 
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A free Jet is usually described in terms of two regions, 
each characterized by different cross sectional velocity profiles 
(ref. 10). The first region, known as the zone of establishment, 
begins at the jet orifice and may be idealized as a jet core, 
.characterized by a flat velocity profile, surrounded by a turbulent 
fluid. The core is rot^hly conical in shape, and diminishes in 
cross sectional area along the jet centerline as a result of shear 
produced by the differing mean velocities of the jet and the surround- 
ing fluid. When the core disappears (at some S/D defined as the 
critical length), the. flow is said to be fully established, and 
is characterized by a velocity profile roughly Gaussian. in nature. 

Albertson (ref. 10) conducted a theoretical analysis of a 
■free jef and formulated relationships for the entrainment in the 
zone of establishment and in the region of established flow. 

Albertson found that the ratio of cross sectional jet volume 
flux to volume flux at the jet orifice increased in a parabolic 
form 

= I to.o63(§)-i-o.oia8(f) . ( 12 ) 

Qo 

in the zone of establishment. In the zone of established flow, 
the relationship is linear 

t = 0.32(f) 

It is convenient to express the amoimt of entrainment in terms of 
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an entrainment coefficient E 

C - X JQ 
^ ~ 06 (KsId) 

Using Albertson’s relationships, the entrainment 
be expressed • 

0.083 -i- 0,o^ 56 

in the zone of establishment, and 

£ ^ 0.3Z 

for established flow. The constant entrainment coefficient 
of 0.32 for a fully developed jet has been confirmed by Ricou 
and Spaiilding (ref. 11) and Saha (ref. 15). 

A free jet entrains surrounding flxxid primarily though 
turbulent shear resulting from the difference in velocity of 
the jet and the surrounding fluid. The entrainment- mechanism 
is considerably more conpiex in the case of a jet in a- crossflow. ■ 
Keffer (ref. 14) states that a jet in a crossflow entrains surrounding 
fluid not only through turbulent shear, but also through the effects 
of free stream and vortex upwash components perpendicular to the 
jet trajectory, in addition, increased '’shear at the boimdary of ' 
the jet resulting from the presence of the crossflow results in 
a more rapid degradation of the jet core and a decrease in the 
critical' length. 

Because of the lack of experimental data concerning entrainment, 
by a jet in a crossflow, it is necessary to attempt to qualitatively 


( 14 ) 

coefficient may 


(15) 


(16) 
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estimate the entrainment coefficient E for use in the conqputer 
model. It is logical to assume that the entrainment coefficient 
E will be larger than 0.32, the value for a free Jet. It is assumed 
that E is constant for a fully established deflected Jet, based 
on observations of free Jets. An upper limit for the entrainment 
coefficient has been suggested by Fearn (ref. 13) in an atteng>t 

N. 

to account for the observed jet trajectory entirely on the basis 
of mass entrainment by the Jet. For a velocity ratio of 8, an 
entrainment coefficient of approximately 1.2 was calculated for 
the region of established flow. Therefore, the value of E used 
in the model would be expected to have a value between 0.32 and 
1.2 . 

The strength density of the line sinks is expressed in terms 
of the change in Jet volume per unit length S/D, or 

vik» 

Relations for q for a Jet in a crossflow- were found by assuming 

= (18) 

for the zone of establishment, and 


for established flow. These relations are of the same form as 
Albertson’s equations for free Jets. By definition 


r _ -L cfO 
^ Qo 



( 20 ) 



16 


The coefficients K_ , K_, and K are found by applying boundary 

1 A 

conditions. It is assumed that at S/D=0, =0 at 

S/B=0, ^ "d^D) continuous at the critical length. Equations 

(18) and (19) then become 


Qo ^ ^ 


^ (s^/o) 

for the zone of establishment, and 


(s/D) 


( 21 ) 





( 22 ) 


for established flow. From these relationships, tjie strength 
density q may be witten 


„ EQq(s/d) 
® (SJD) 


in the zone of establishment, and 



for established flow. The amount of entrainment is therefore 
defined by establishing values for E and S^/d. From Keffer (ref. 14) 
and Pearn (ref . 13) , the critical length S^/D is equal to about 3 
for a velocity ratio of 8. This value, although approximate, has 
more experimental basis than values of E for deflected jets. 
Therefore, E was the only parameter which was varied in the model. 
The value for E was adjusted until the model provided good overall 
correlation with the experimentally determined pressure distribution. 

^i'th the preceding description of the vortices and sinhs 
incorporated into the model, it was possible to calculate the 
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velocity induced on the flat plate. This is accomplished through 
vector summation of the velocities induced on the flat plate by 
each incremental vortex and sink filament. image vortex and sink 
distributions are instituted "below the flat plate in order to 
establish the boundary condition at the plane of the flat plate, as 
shown in Figxjre 3, With the velocity field on the flat plate 
determined, the pressure coefficients at each field point may be 
calculated from Bernoulli ’ s equation. The pressure coefficient 
is derived from Bernoulli's equation and may be written 

where U is the fluid speed at a particular field point. 

The pressure distribution on the flat plate is most conveniently 
presented by defining the field points in terms of polar coordinates. 
Figure 4 shows the tunnel and field point coordinate systems used in 
the experimental investigations of Fearn and Weston (ref. 3). Each 
field point is expressed in terms of r/D and 9. Fearn and Weston 
presented pressure coefficient data in terms of contour plots and 
pressure distribution along rays (constant angle 9) and circles 
(constant r/D). The results of Fearn and Weston are used for 
comparison with the thecMcetical results of the model in Figures 5 
through 16. 

As expected, the pressure coefficients predicted by the 
model were found to be in considerable error in the wake regipn, 
due to viscous effects and separation. For the purposes of this 
- study, a wake boundary is defined by con^jaring contour plots of 
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experimental and theoretical pressure distributions and fitting 
a power curve through the points where the theoretical results 
begin to deviate significantly from the exj>erimental contours. 

This wake boundary, which is defined for the purposes of this study 
only, should not be confused with that which would be observed 
through oil smear studies or other flow visualization techniques. 

It is assumed that a point on the wake boimdary defines the 
pressure at all points in the wake region with the same X-coordinates 
as the point on the boundary. The elimination of Y-dependent 
pressures results in straight contour levels extending from the wake 
boundary to the X-axis. From studying contour plots of the eotperi- 
mental pressure distribution for various velocity ratios (ref. 3), 
it appears that this method for approximating the pressures in the 
wake region should be applicable for velocity ratios greater than 5. 
At lower velocity ratios, the contours in the wake region cannot 
be approximated by straight lines, due to the tendency for the 
contours to become increasingly curvilinear at lower velocity ratios. 

Lift and pitching moments on the flat plate are calculated 
for the complete model. Lift is calculated by summing forces on 
incremental areas of the flat plate. Moment about the Y-axis is 
found by multiplying incremental areas by the distance from the 
Y-axis, and summing over a large area of the flat plate. 



RESULTS 


The simplest form of the model consists of the vortex distribution 
alone, with the vortices intersecting the flat plate at the center 
of the Jet orifice. In this simplified form, it is assumed that the 
flat plate pressure distribution results from the vortex pair and 
free stream interaction only. Figures 5 through 8 compare the 
experimental pressure coefficients with the theoretical values 
predicted for this case. Figures 5 and 6 show the variation of the 
pressure coefficient with radial distance from the center of the Jet 
orifice for given values of 6. The plots show fairly good agree- 
ment with experiment from about 9=75° to 0=120°. From 120 to 180 
the theoretical values deviate considerably form the experimental 
results, as expected, since conditions in the wake region render 
potential flow considerations invalid. The theoretical-values from 
0=0° to 0=75° also deviate significantly from the experimental 
values, due to the neglection of Jet flow entrainment, as will be 
shown later. Figxire 8 shows a contoixr plot of the pressure 
distribution induced by the vortices alone. The theoretical 
contours clearly deviate greatly from experimental values outside 
the 0=75° to 0=120° range. 

Figures 9 through 12 show the effects of adding a sink 
distribution along the Jet centerline to account for entrainment 
effects. In this study, good results were obtained for E=0.6, 
with a critical length of 3 Jet diameters. The value of the 
entrainment coefficient is within the range expected, i.e. 0.32<E'^L.2 . 
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It is evident that the addition of a sink distribution greatly 
improves the accuracy of the model in regions upstream of the Jet 
orifice. The greatest errors outside the wake region occur at ■ 
low r/D (<1.0) during the transition from positive to negative 
pressure coefficients in the range 0=20° to 0=45°, as can be seen 
in- Figure 9. However, Figure 12 shows that the theoretical and 
experimental contours are quite close in this region. It may also 
be -seen that the density of contours in this region indicate .that 
high pressure gradients exist. Hence, a small displacement of the 
contours results in a large change in pressure values displayed in 
the radial plots, suggesting that' the rather large. errors observed 
in the radial plots for this region are a result of the manner 
in which the data is presented, rather than a result of a serious 
failiire of the model; Good agreement is obtained up to about 
0=135°, at which point viscous effects and separation appear to 
become important. The largest error observed on the contour 
plot (Figure 12) is on the order of one Jet diameter, and involves 
the -0.2 contour at about 45°. This is not considered to be a 
particularly significant error, since insignificant variations in 

■ i 

pressure can cause noticeable shifts in the placement of contours 
in regions of very small pressure gradients. 

The potential flow model results in pressure coefficients 
being too large in the wake region. As stated earlier, the failure 
of the model in the wake region is not unexpected. However, a 
complete description of the pressure distribution must take the 
wake region into account in some manner. Figures 13 through 16 
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show the resulting pressure field for a model incorporating both 

jet entrainment and wake region considerations described in* the 

previous section. It may be seen that the accuracy of the 

theoretical wake region pressures is greatly improved, althoTi^h 

good quantitative agreement has not been* acheived for r/D of less 

than 1.5 . Figure 16 shows the contour plot for the complete 

model. From the contotir plot, it may be seen that the largest 

discrepancies occur in the wake region for C =-0.1 . This is not' 

P 

considered to be significant, since this region is characterized 

\ 

by very low pressure gradients and a large uncertainty in contour 
placement. Figures 17 and 18, show the theoretical results of the 
complete model compared with the experimental results of Thontpson 
(ref. 4) and Bradbury and Wood (ref. 5). It may be seen that the 
contours agree quite well with Thompson's results. Agreement with 
the results of Bradbury and Wood is not as close, but is still 
considered to be within the range of experimental error. 

- The lift and pitching moments on the flat plate are of special 

f 

interest, particularly in the case of transitioning- TOOL craft. 
Figure 19 shows a plot of experimental values of M/^D and I/T 
versus the velocity ratio R, The circles represent values of 
M/TD and 1 /t predicted by the complete model for a velocity ratio 
of 8. The theoretical values show good agreement with the 
experimental values. 

Sellers (ref. 7) formulated a model in which' the vortices 
intersect the flat plate at the edges of the jet orifice, rather 
than at the center. Although this model may seem more physically 
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reasonable, the resulting pressure distribution predicted by the 
model (Figure 20) deviates considerably from the experimental 
data. The .extremely low pressiire coefficients result from a 
combination of’ the effects of non-zero vortex strength at the 
plane of the flat plate and the proximity of the filament vortices 
to the field points closest to the jet orifice. It is evident that 
a model based on vortices beginning at the edges of the jet orifice 
is incompatible with the use of filament vortices. 



CONCLUSION 


The pressure distribution calculated from the model presented 
in this thesis is in fairly good agreement with oxporimental data. 
The use* of this model was demonstrated for a velocity ratio of 
8 in order to make use of an available analytic description of the 
characteristics of the vortex pair. The model is equally applicable 
for other velocity ratios when used in conjunction with graphical 
descriptions of the vortex properties. The method for predicting 
wake region pressures is expected to be applicable for velocity 
ratios greater than 5; at lower velocity ratios, curved lines may 
be necessary to approximate the shape of the wake contours. 
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Figure 1. Sketch of VTOL Aircraft Transitioning 
from Hovering to Horizontal Flight 
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PiguSre 3. Arrangement of Singularity Distributions 

in the Model 
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Figure 4. Tunnel and Field Point Coordinate Systems 



Figure 5, Comparison of Model Results with Experimental Pressure 
Coefficients for Zero Entrainment (0 <0<9O ) 
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Figtire 6. Compaxison of Model Results with Experimental gressure 
Coefficients for Zero Entrainment (105 <0<18O ) 
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■Figiire 10. Comparison of .Model Results with Experimental Pressiare 

Coefficients for E=0.6 (105*^9<180°) 
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•Figvire 11. Comparison of. Model Results with Experimental Pressure 

Coefficients for E^.6 ( .58 -^r/D 1.00) 





Figirre 12. Compaxison of Model Results with Experimental Pressure 

Contours for E=0.6 
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o RIGHT HfiLF-PLflNE (9) 
^ LEFT HfiLF-PLfiNE (-9T 

model 



Figure 15. comp^ison of-Mod&l Results with Experimental Pressi^e 
Coefficients with Application of Wake Region 
Considerations (E=0,6, ,58 < r/D ^ 1.00) 
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Figure 16. Comparison of Model Results with Experimental Pressvire 

ContovjTs with Application of Wake Region Considerations 

(E=0.6) 
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Figure 18. Comparison 


of Model Results with Results of Bradbiiry and Wood (ref . 5) 
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Figure 19. Comparison of Model Results and Experimental Lift and 

Moment on the Flat Plate 




Figure 20. Comparison of Model Results and Experimental Pressure 
Coefficients (Vortices' Begin at Edges .of Jet Orifice) 




APPENDICES 



APPENDIX I 


Derivation of Finite Filament Sink (Source) Flow 


The computer model institutes a distribution of finite 
filament sinks along the Jet centerline. It is necessary to 
•formulate an expression relating the sink filament strength 
and length to the induced velocity at a field point, consider 
a distribution of equal strength sinks along the X-axis of a 

i - X 

■ local coordinate system. (The axes used in this derivation 

should not be confused with the tunnel coordinate system described 


in- the main text.) 



The total velocity potential induced by a line sink of length b 
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may be written 


4 > = 



The X-velocity induced may be foimd by the expression 

.t 






Therefore 


Let 








x-i> 




‘"^Jcr ^ yT‘ 

X 
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^ y 

X-i) 

Vtry L [(X‘i>)\f y*J f X* y*J 


-I —2— r Coi" o<a " COS <^2 7 

^ Hirh L J 


The problem essentially reduces to a two dimensional 
geometrical analysis. When the endpoints of the filament and 
the field point are defined, the values of AP, OP, h, and 
may be found through planar geometry. The projections of 

2 

U and U onto the flat plate will yield the velocity, on the 
X y 

flat plate which is induced by the finite filament. 



APPENDIX II 


Derivation of Finite Filament Vortes Plow 


The computer model institutes a distribution of finite 
filament vortices along the vortex trajectories. Although the 
vortex endpoints and the field points are specified in three 
dimensions, the problem is two dimensional as far as calculating 
the induced velocities on the flat plate is concerned. Consider 
a filament 'Cortex aligned along the X-axis of a local coordinate 
system (not to be confused with, the tunnel oriented coordinate 
system) , The X-Y plane contains the field point and the two 
vortex endpoints. 



The general equation for the velocity induced by a filament vortex is 



where P Is the strength of the vortex. In vector form the terms 
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in the general equation become 


5 = (X,-X')t + y.j 

dI(x') = oli't 
S X d/(xj = y. dx' k 


Let 


^ = Xo-X' 


Therefore 




U = 


Jl 

Hit 


n. 


X„-Xa 




JL f 




[U 


X^-X, _ ] 


- J2^ fcas - cos &g] k 
t^Trh ^ 


If the. field point and the endpoints of the vortex filament 
are defined, the parameters h, 0 , and 0„ may be found throiigh the 

\ Jm ^ 

use of planar geometry. The projection of U onto the plane of the 
flat plate will equal the velocity induced on the flat plate by 
the filament vortex segment. 
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APPENDIX III 


COMPUTER PROGRAM 

FOR CALCULATING PRESSURE DISTRIBUTION 


ON FLAT PLATE 
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c 

C THIS PROGRAM CALCULATES THE PRESSURE COEFFICIENTS ON'THE FLAT 

C PLATE INDUCED BY A JET ISSUING NORMALLY INTO A CROSSFLOW* 

c 

DIMENSION X{200),Y('200>,Z(2CO)tCP(200), SI200) 

DIMENSION RADC26J ,THETA(15) 

DIMENSION P126,15) 

CaMMaN/ONE/GAM,CONPT,CONST 
COKMCN/TWO/RAO, THETA rP 
COMMON/THREE/X, Zt S»N 

COMMON/FGUR/XO,YOTZO,XXtYltZi,X2,Y2,Z2t VXTQT,VYTCT» VZTOT 
COMMON/F IA/E/N0FIL»R7VJX,VJY,VJZ 
INTEGER OPTION 
COMMON/SIX/OPTION 

DATA RAD/0.58,0.62,0.66,0.70»0.75T0.80,0.90Tl*00»1.20,l.AOf 1.60» 
11.8Of2-OO»2.25,2.5O,2.75»3.G0» 3. 25, 3,50,3.75,4.00,4.50,5.00,5.50, 
26.00,7.00/ 

DATA THETA/0.0, 10.0,20.0,30.0,45.0,60.0,75.0,90.0,105.0,120.0, 

• 1135.0# 150. 0,160. 0,1 70. 0,180.0/ 

AE=00-3473 
BE= 1.127 » 

C^=6.4291 

DlA=0i.333-33 

U=127.0 

READ 15,66) OPTION 
C . 

C .• IF QPTION=l, THE VORTICES ■ START FROM THE CENTER. OF THE JET ORIFICE. 

C IF aPTI.QN^2, THE VORTICES START FROM THE EDGES OF THE JET ORIFICE. 

C , ■ - . 

READ 15,65) NOFIl 
READI5,24) R 
NOWAKE=l 
C 

C IF NCWAKE EQUALS 1, SUBROUTINE WAKE IS CALLED. 



c 

N0SINK=1 

C IF NOSINK EQUALS 1, A SYSTEM OF LINE SINKS IS GENERATED ALONG THE JET 

C CENTERLINB. QTHERWISEt NO JET ENTRAINMENT IS ACCOUNTED FOR* 

C 

6AM=0*0 

DO I ICO M=;lt-15 
DO 1100 K=l,‘26 
P.(K,M)=1E6 
1100 CONTINUE 
N=N0F1L 
CALL ENDPT 
CPUJfcO.O 

IF CORTICN.EQ.n CONST=2*OA 
IF (0P.TICN.EQ.2) C0NST=1*389 
00 17 

00 17 K=l,26 
ARC=0-0 

THET~’THETA(M L*3. 1416/180*0 
C 

C DEFINE FIELD POINTS. 

C 

XO=-RAO(K)*COS(THET) 

YO=RAO<K)»SIN(THET> 

' Z0=0.0 
VXT0=0.0 
V.YT0=0.0 
VZT0=0.0 

IF(OPTION.EQ.l) Y(il)=;0.0 
IF(0PTI0N’.EQ.'2) YU)=i0.5 

IF (N0WAKE.EQ.I.AN0.THETA1M).GT.90) CALL ViAKE(K,M) 

DO LOOP 1 GENERATES A SERIES OF VORTEX FILAMENT SEGMENTS IN ORDER 


C 

c 
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c TO APPROXIMATE A CURVED VORTEX FILAMENT. 

C 

DO 1 1=1, N 
X1=X(I) 

X2=X(X+U 

SEG=SU> 

ARC=ARC+S( I) 

cpn+i)=ARC+sa+i)/2<:o 

CONPT=CPU) 

H0=CGNST*ll.-»EXPl-ARC/8.0) ) 

IF ( OPTION. EQ.:2) H0=H0+0.5 
BD=2.11/S0RT(ARC) 

ETA=BD*H0 

YI I+1)=H0/ERFIETAT 
Y1=Y(I) 

Z1=ZII) 

Y2=Y<I+1) 

Z2=Z{I+1) 

C 

C DC LOOP 4 GENERATES TWO OBSERVED VORTICES AND TWO IMAGE VORTICES.' 

C ; 

DO 4 J=l,4 

If (ZO.'EQ*0.0.:AND.J.GE.3) go TO 4 
GO TO C9,6,7»6),J 

6 Y1=-Y1 
Y2=-Y2 
GO TO 9 

7 Z1=-Z1 
Z2=-Z2 

9 CONTINUE 

c ■ 

C CALCULATE DISTANCES FROM VORTEX SEGMENT ENDPOINTS TO FIELD POINTS. 

C 

A=SQRTnXl“XOl**2+(Yl-YO)**2+(Zl-ZO)=«‘=»=2) 



• B=nSQRT( (X2-X0)=**2+( Y2-Y0)**2 + tZ2“Z0)**2 ) 

C=SQRT { ( X2-X1 )'**2+( Y2-Y1 1**2+{ Z2— Zl ) 

T=(A+8+C)/2. 

B«OUT=T*(T-Ay«(T~B)*(T-C> 

IF (BMOUT.LT^O.O) GO TO 502 
E=SQRT(BMOUT) 

GO TO 503 

502 E=^0.01 

503 CONTINUE 

PERPENDICULAR DISTANCE FROM VORTEX FILAMENT TO FIELD POINT. 

H^Z.^E/C 
H=H*BIA ■ 

IF {Hi,LT;.00ir H=-001 

C0S1=( A**2+C<5*2-B**2)/I 2 .0*A*C) 

C0S2=>I A#=P2-8**2-C#*2)/( 2-0*B#C) 

IF {R.EQ--8.0) CALL GAMMA8 
G1=2.0*GAM*U«DIA 
IFU.EQ.Nl CBS2=-1 

FIND MAGMITUDB CF VELOCITY INDUCED AT FIELD POINT BY VORTEX SEGMENT. 
V=GI*ICOS1-COS2)/(4.0*3. 1416=>H) 

CALCULATE VECTOR COMPONENTS FROM CROSS PRODUCT OF VECTOR A AND VECTOR B 

AC= ( Zl-ZO 1=*^ (.Y2-Y0 )- 1 Z2-Z0 1’J' ( Y 1- YO ) 

- 8C=(Xl-XO)*(Z2-ZO)-(X2-XO)*(Zjl-ZO) 

CC= ( Vi-YO )*( X2-X0 ) - ( Xl-XO J * ( Y|2-Y0 ) 

GO TO (14,15,14,15) ,J 
15 CONTINUE 
AC=-AC 
BC=-BC 
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14 CONTINUE 

C FIND HIRECTIGN COSINES OF INDUCED VELOCITY VECTOR. 

C 

D=SQRT { AC*<‘2-^8C«*2+CC=«‘*2 ) 

DC0SA=AC/0 

DCOSfi=BC/D 

DCOSG=CC/D 

C 

C SUM INDUCED VELOCITIES- 

C 

V XT 0= V XT 0 +D COS A =» V 
VYTO=VYTD+DCOSB=<‘V 
VZTO=-VZTC.+OCOSC*V 
4 CONTINUE 
1 CONTINUE 

IF{ZO-.EQ.O.O) VXTO=2.0*VXTO 
IFCZ0i.EQ*O.O) VYTO=2-0’«'VYTO 
if(zc.:eq-o-o) vztg=o.o 
VJX=0-0 ' 

VJY=0-0 

VJZ=0.0 

IF (NOSINK.EQ.l) CALL JET 

C ADD VELOCITIES INDUCED BY SINK FILAMENT SEGMENTS. 

C 

VXTO=VXTO+VJX 

VYTO=VYTO+VJY. 

VZTO=VZTO^tVJZ 

VXTO=VXTC+U 

I FIND TOTAI. VBUQCITY AND PRESSURE COEFFICIENT ON THE FLAT PLATE 

c 



noon 


VTOT=SQRT (VXT0*4t2+VYT0**2+VZT0**2) 

P<Kt«)=l*0-( VTOT«*2.0)/(U**2.0) 

17 CONTINUE 
CALL 6XPCP 
CALL PRTMAT 
CALL LIFT 
2A FQRHAT(F5.2> 

65 FORMAT (13) 

66 FORMAT (ID 
STOP 

END 

SUBROUTINE ENDPT 

SUBROUTINE ENDPT GENERATES ENDPOINTS IN X AND Z FOR EACH VORTEX 
SEGMENT. 

COMHON/THREE/Xt 200)tZ(200)tS(200) ,N 

AE=0.3473 
BE=1.127 ■ 

CE=0.429l . 

R^8 • 0 t 

Al=CE*AE*{R«*BE) 

zm=o.o 

Z(2)=.025 

X(1)=(Z( 1)/(AE*(R**BE)) )*«( l.O/CE) 

X(2)=(Z(21/( A6*IR**BE)) )*#( 1.0/CE) 

S(1) = ((X(2) — X(1))**2.0+(Z(2)— Z(l) )**0. 5 

SCNE=2.0*S(U 

NN=N+1 

T0TARC=0.0 

DO 100 I=2«NN 

02=Al*tX( n*=«KCE-l.on 

PHI=ATAN(DZ) 

X( I + D=SONE4'COS(PHD -hXI I ) 
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Z( I + 1)=AEj^(R*#BE)*(X( I + l) )**CE 
CALCULATE LOCAL RADIUS OF CURVATURE. 

RADCUR=( l.:08375’M XII ),**!. 047267 )^tl.+2. 41041=5' (Xt I ) ** |-1 . 1418 ) ) ). ) *^ 
11.5 

. Sm=SQRT(|X( I + 1)-*X( I) >4=«‘2 + (Za+l)-2( I) )=4'=4'2) 

TQTARC=TCTARC-»:Sn ) 

CALCULATE LENGTH OF NEXT FILAMENT BASED ON RADIUS OF CURVATURE. 

IFCTQ/TARC.GT.l.OO) SGNE=RADCUR/50. 

WRITEC6,101) TOTARCiRADCURfSONE 
101 F0RMAT{5X,«T0TARC,.RADCUR,S0NE ='t3F10-4) 

100 CONTINUE 
RETURN 
END 

SUBROUTINE JET ■ 

SUBROUTINE JET INSTITUTES A SERIES OF LINE SINKS ALONG THE JET 
CENTERLINE IN ORDER TO ACCOUNT FOR MASS ENTRAINMENT BY THE JET. 

DIMENSION X<200)#ZI200),CP(200)fS<200> 

COMMON/FOUR/XO,YO,ZO*AZ.BZ,CZ*ZZ,EZ,F2,HARPO,ZEPPQ, GROUCH 

COHMCN/FIVE/NOFlL»RtVJX,VJY»VJZ 

COMMON/SEVEN/E, SCR IT 

cpm=o.o 

AE=0.9772 

BE=0.9113 

CE=0-3346 

U=127.0 

DIA==.3333 

PI=3. 14159 
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C 

C 

C 


C 

c 

c 


c 

c 

c 

c 


*E* IS THB ENTRAINMENT COEFFICIENT. *SCRIT» IS THE CRITICAL LENGTH FOR 
ESTABLISHMENT OF THE JET FLOW. 


6 =. 0.6 

SCRiT5=3.0 . 

GENERATE THE ENDPOINTS OF THE SINK FILAMENT SEGMENTS. 


XT0T=O.0 


ZT0T=0.0 


T0TARC=0.0 

Z(1)=0.0 


Z{2) = .‘025 

XU ) = IZI 1.K/I AE^tR*=«‘8E )))*♦( 1.0/CE) 

X(2> = (Z(2I/< AB«(R**8En )♦*( i;0/CE) 

SI'l ) = t(X ) )**2.0+ (Z(2>-2( 1 ) )**2.0 )#=«=0- 

XT0T=Xt2)-X( 1) 


5 


SONE=2-0*SI U 


NN=NCFIL+T 
ARCLTH=0.O 
DO lOG 1^2, NN 

DZ=CE«AE^IR#t*BE)*<XTOT**(CE-1.0n 

PHI=ATAN(DZ) 

X( I + l)=SONE*COS(PHl HiXTOT 
z(i+i)=AE+iR«i»BE)=«‘(xn+i) )=^*ce 
S(I)=SQRTHXH + X)-XI I ) ) =>’(‘2+1 Z (I + 1 )“Z U ) )**2> 
XTOT=xnU) 

ARCLTH=S (I) +ARCLTH 


CALCULATE LOCAL RADIUS OF CURVATURE OF JET TRAJECTORY AND ADJUST LENGTH 
OF NEXT SEGMENT ACCORDINGLY. 


IF(XM).LT.O.’Ol) GO TO ICO 

RA0CUR=i.4475=MX( I)**(~1.66 54) )/U.+4.73l5*<X(I )=>=> 


(-1.3308) )T**1.5 


C7I 

00 
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RA0CUR=1.‘/RADCUR 

IF(ARCLTH«GT*3.00) S0NE=RADCUR/50. 

100 CONTINUE 
VJX=0«0 
VJY=0-0 
VJ2=0.0 

DO 13 I=L,NOFIL 
TOTARC=TOTARC+S(I ) 

CPC I+1)=TQTARC+S( 1+11/2.0 
CONPT=XP( l> 

Xl=X(U 
X2=X C I+l ) 

Y1=0.0 

Y2=0-0 

zi=zm 

Z2=ZU+1) 

DO 19 J=l’,2 

IFCZO.EQZO.O.AND.J.EQ.Z) GO TO 19 

IF (J.EQ.l) 00 TO 21 

Z1=-Z1 

I 

21 CONTINUE 

CALCULATE STRENGTH OF SINK FILAMENT SEGMENT. 

UJET=R*U 

QO=UJET=fr(RI«iOIA*#2)/A.O ) 

IF (CONPT.LEZSCRIT) Q==E*eONPT*QO/ ( SCR IT’9=0IA ) 

IF (CONPT\GT.!SCRIT) Q=E*QO/DIA 

CALCULATE THE DISTANCES FROM THE ENDPOINTS TO THE FIELD POINT AND 
THE HEIGHT OF THE RESULTING TRIANGLE. 

PA=SQRT C C XO-Xl ) ♦^Z+C YO-Y 1 >=*'=9=2+ 1 ZO-Z 1 ) **2 ) 
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PB=SGRT{ (XO-X2)*=«‘2 + (y6~Y2)*=4'2 + ( 10-I2)^^2) 

AB=SGRT( (Xl-X2)**2+( Yl-Y2)«*2+( Zl-Z2)^*2 ) 

SL=(PA+PB-+AB)72.0 

ARZ=SL«( SL-PA) * ( SL-P B )=*' ( SL-AB ) 

IF (ARZ,LE.O.O) GO TO 301 
AK=SQRT(ARZ) 

QP=2..0«AK/AB 
GO TO 302 

301 OP=.Ol 

302 CONTINUE 

AQ=SQRH (PA*«2.0W(QP**2.0} ) 

C0Sl=(PA*=*2+A.B*=»2-PBSi#2 )/(2.0#AB’«'PA3 
CQS2 = ( AB*^2 + PB**2^PA=<5=«‘2 )/{ 2-0*AB=«'PB ) 

IFICQSl.GfT.l.O) COS 1=1.0 
If <C051.LT.-1.0 ) C0S1=-1.0 
IF(C0S2.G>T.1.0) COS2=1.0 
IFICCS2.LT.-1.0) C0S2=-1.0 
C0S2=~C0S2 
ANGI=ARCOS(COS1 ) 

ANG2=.^RC0S<CCS2> 

C . '* 

C CALCULATE THE DIRECTION COSINES OF THE FILAMENT SEGMENT. 

C 

CQSALF=(X2-XU/AB 
C0SBET=IY2-YH./A8 
C0SGAM=(Z2-ZL)/AB , 

C CALCULATE THE X» Yt AND Z CO-ORDINATES OF .THE POINT . OF INTERSECTION 

C Of THE PERPENDICULAR LINE FROM THE FIELD POINT TO THE FILAMENT. 

C 

IF (ANGl.LT.Pl/2.) GO TO ll 
IF IANGltfGT.PI/2.) GO TO 10 
11 XQ=X1+AQ#C0SALF 
YQ=YH-AQ*COSBET 
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ZQ=Z1+AQ*C0S6AM 
GO TO 12 

' 10 XQ=X1-AQ’^C0SALF 

YQ=Yl-AQ=«lCOSBET 
Z.Q=Z1-AQ^C0SGAM 
12 CONTINUE 

CALCULATE THE DIRECTION COSINES OF THE PERPENDICULAR LINE FROM 
THE FILAMENT TO THE FIELD POINT* 

COSPH'I = { XO-XJQ)/GP 
C0S0MG=( Y0-YQV./CP 
COSZET=(ZO-ZQ)/QP 
QP=QPs«‘DIA 

CALCULATE THE VELOCITIES INDUCED BY THE FILAMENT SEGMENT AT THE 
FIELD POINT. 

UY=Q*<C0S2-C0S1 }/ (4.0*P I*QP ) 
ux=Q«ni.a/PA»-a.o/PB) )/(a.o*pd 

UX=UX/DIA I 
VXSINK=UY*COSPHI 
VYSINK=UY=0fCOSQH6 
VZSINK=UV<!COSZ'ET 
UXSINK=UX.*COSALF 
UYSINK=UX»COSBET 
U2SINK=UX^C0SGAK 

ADO THE RESULTING VELOCITY COMPONENTS^ 


o> 

M 


VJX=VXSINK+UXSINK+VJX 
VJY=VYSINN+UYSINK+VJY 
VJZ=VZSINK+UZSINK+VJZ 
19 CONTINUE 



c 

c 

c 

c 


c 

c 

c 

c 
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CONTINUE 

IF ( ZO.EQ.0.0) VJX=2.0«VJX 
IF(ZO.EQ-O.O) VJY-2.0*VJY 

if(zo«eq.:g-o) vjZ"0-o 

RETURN 

END 

SUBROUTINE GAMMA8 


SUBROUTINE GAMMAS SELECTS AN 
FOR EACH VORTEX FILAMENT FOR 


appropriate vortex strength 
A VELOCITY RATIO OF 8.0 . 


COMMCN/ONE/GAM,CONPT, CONST 

INTEGER OPTION 
C0MMGN/SIX/0.PT10N 

•CCNPT* IS THE DISTANCE ALONG THE VORTEX TRAJECTORY TO THE CONTROL POINT 
OF THE VORTEX SEGMENT BEING ANALYZED. 


IF ICQNPT.'EQ.O.O) GO TO 103 

BD=2-ai0268/SQRT(C0NPT) 
HO=CONST*U.-EXPI-CONPT/8^0) ) 
IFtOPTION-EQ.2) H0=H0*0.5 
ETA=BD#HO 

GAM=0.72*-8*0*ERFt ETA) 

GO TO 104 

103 IF {OPTION.EQ.'l ) GAM=?0-G 
1F(CPTI0N-EQ*2) GAH~5«76 

104 CONTINUE 
RETURN 
END 

SUBROUTINE PRTMAT 


C 

C 

C 


SUBROUTINE PRTMAT PRINTS CALCULATED COEFFICIEN^^ 

EXPERIMENTAL PRESSURE COEFFICIENTS* THE DIFFERENCE BETWEEN 


01 

to 






C CALCULATED AND EXPERIMENTAL PRESSURE COEFFIC IENTS» AND UNWEIGHTED 

C STANDARD DEVIATION- 
C 

DIMENSION DIFFI26t15I 

COMMON/TWG/RADC 26 ) f THET A ( 15 ) tP( 26« 15) 

COMMCN/FIVE/NOFIL^R»HUPfTUP»THREEP 

INTEGER OPTION 

COMHON/SI'X/OPTION 

COMMGN/SEVEN/E,SCRIT 

CQMM0N/TEN/D(21,26) 

DELPT=0-0 
SUMDIF=0.0 
DO 999 M = l,1.5 
DO 999 K=?l,26 
C 

C CALCULATE. DIFFERENCES IN EXPERIMENTAL AND THEORETICAL VALUES. 

C ■ 

. DIFF{K,M)=P(K,M)-D(M^K) 

IF (ABS(DIFF(K,M)).GT.10) D I FF ( K,M ) =1E6 
' IF(DIFF( KrH) ^EQ.1E6) DELPT=D£LPT+1 . 

IFIDIFFIK^M) .JEQ.IE6) GO TO 999 
DIFF{K,Mls=DIFF< K,M)=i'.*2 
SUMDIF=SUMDIF-»:DIFFI KtM) 

C 

C CALCULATE. PERCENT DIFFERENCE FOR EACH POINT. 

C 

DIFF{K,M)=SQRT(DIFF{K^M) )*1C0. 

999 CONTINUE 
C • 

C CALCULATE UNWE IGHTED' ST ANDARO DEVIATION. 

C 

STDEV=SQRT{ { SUMDIF/ ( 390.-DELPT) ) ) 

WRITE(6,101l); NOFILiR 
WRITE(6,1013) EtSCRIT 


o 

CO 



1004 


1005 

1000 

1001 

1002 

1003 

1006 
1007 

1009 

1010 
1014 


IF (CPTIGN.EQ.I ) WRITE(6,1009> 

IF {OPTION. EQ. 2 ) MR ITE( 6 , 10 10 ) 

WRITE (6tl012). 

WRIT£{6,1007) 

WRITE(6,1001) 

ViRITE{6,1002) {RA0{K),,.K=*1, 13> 

WRITE<6, lOOir 
WRITE(6» 1006). 

DO 1004 M=l,15 

WRITE (6, 1003) THETA(M)» ( P ( K , K=1 , 13 ) 
WRITE(6, 1014) ( DiM,LltL=l,13> 

WRITE! 6*1015). (CIFF{K*M)*K=1*13) 
CONTINUE 
WRITE{6, 1000) 

WRIT£<6* 1007) 

WR1TEC6,1001) 

WRITE<6,1002) 

WRlTe{6,1001) 

WRITE{6, 1006) 

DO 1005 M=il*15 

WRITEI6, 1003) THETA! M), ( P { K ,M ) , K=14,26) 
(D!W,Ut.L = 14,26) 
(0IFF(K,M)tK=14,26) 


!RAD(K),K=14*26) 


COEFFICIENTS ,///) 


WRITET6, 10141 
WRITE<6, 1015) 

CONTINUE 

WRITE(6, 10161 STDEV 
F0RMAT(1H1,///»55X,21HPRESSURE 
FORMAT! 1X»-127{1H* ) ) 
F0RMAT(15X,13(F5.2,4Xn 

F0RMAT!3X,F4.0,2X,1H*».5X,13! F7.'3*2X) ) 

F0RMAT(4X»*ANGLE*,/) 

F0RMAT!64X* »R/D* ) 

FQRMATC5Xt*0PTI0N=lA, VORTICES START FROM 
FORMAT(5X»’OPTION=2A, VORTICES START FROM 
FORMAT (1X*» EXP. VALUES* »3X» 13( F7.3j 2X) ) 


CENTER OF JET ORIFICE.* 
EDGE OF JET ORIFICE.* ) 



o O O O o o o 


1015 

1016 
1011 

1012 

1013 


FORMAT ( IX t 'PCT« ERROft* »4X»13(F7.3t2X) »/) 

FORMAT t IHl* 5X»* STANDARD DEVIATION =’»F10.4»/) 

FORMAT ( IHlvSXf I3».1X, ‘VORTEX AND SINK FILAMENT SEGMENTS ARE USbC’t 

I/, 5X, ‘VELOCITY RATIO =»,F5.2fy/) 

FORMAT {///,.55X,* PRESSURE COEFFICIENTS*,///) .-ott-ttai j fncth 
FORMAT 15X>* ENTRAINMENT COEFFICIENT =' ,F6.'=i, /,5X, CRITICAL LENGTH 

1 FOR JET FLOW ESTABLISHMENT , F6*4, IX, • S/D • ) 

RETURN 

END 

SUBROUTINE LIFT 

SUBROUTINE LIFT CALCULATES LIFT AND PITCHING MOMENTS ON THE FLAT 
PLATE FROM BOTH EXPERIMENTAL AND CALCULATED PRESSURE COEFFICIENTS. 


DIMENSION AZIMTHI30) 

COMMGN/TWO/RADt 26 ) ,THETA ( 15 ) , P ( 26,15) 
COMMCN/TEN/O(21,26) 

RA0IAN=3. 14159265/ 180.0 
QINF=19-935 
VJ=0-1247 . 

WP=6.66 

WMACH5=0.9-3 

T=148.’10 

TE=(T+459-7)/I1.0+0.2*WMACH«*2) 

CALCULATE THRUST OF JET. 

THRUST=WP-/32.17=^WMACH*49J02=«'SQRT(TE) 

EXLIFT=0-0 

EXPPCH=0*C 

OLIFT=0.0 

PM0M=0.0 

R=4.0 

TAREA^O.O 


C5 

tn 



209 


C 

C 

.C 

C 

C 

C 


c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


00 -209 Mrl*15 

AZIKTH(M)=iTHETA(M):*RADI AN 

00 210 K=l,23 
DO 210 M=I,1A 

DTHETA=AZIMTM{M + 1)'-AZIMTH1M) 

CALCULATE INCREMENTAL AREA. 

AREA=R^R*!(RADIK + 1)«*2-RAD(K )*DTHETA/2. 

CALCULATE AVERAGE PRESSURE COEFFICIENTS OVER AREA. 

CPAVG=(P (K*M) '♦jPtK+l tM).+ P <K»M + 1 ) +P(K + 1»M + 1 n /4.0 
0ELFQR=CPAVG^QINF=»AREA/1A4.0 

CALCULATE MOMENT ARM.. 

PMARM=R=«‘<RA0(K+1 ) +RADIK ) )*CQS ( AZ IMTH ( M+ 1 ) -DTHE TA /2 , ) /2. 
PMDM1=PMARM=«'DELF0R/12. 

SUN LIFT OVER AREA. 

DL I FT=OL IFT+DEL FOR 
SUN MOMENTS. 

PM0M=PM0M,+PM0M1 

SUM AREA. NOTE TOTAL AREA INCLUDES ONE-HALF OF FLAT PLATE. 
TAREA^TAREA+AREA 

CALCULATE AVERAGE EXPERIMENTAL PRESSURE COEFFICIENTS AND SUM LIFT AND 
MOMENTS. 


05 

o 



& 


EXPAVG-(D;(M♦K) + D(M+l^K)+D(MfK+l)+D{M+l•»K+l) )/4«0 
EXPF0R=EXPAV6*Q INF^AREA/ 144.0 
EXP«Cf‘=PHARM*EXPFQR/l2. 

EXLIFT=eXLIFT+EXPFOR 

EXPPCH=EXPPCH+EXPMOM 

210 CONTINUE 

C BY SYMMETRY, CALCULATE TOTAL LIFT, MOMENT, AND AREA BY MULTIPLYING ABOVE 

C RESULTS BY TWO. 

C 

DL IFT’=2.0=»‘DL IFT 
PMGM=2.0*PMOM 
TAREA=TAR£A*2.0 
EXLIFX=2.0=5'EXLIFT 
EXPPCH=2.0#6XPPCH 
C 

C NON-DIMENSIONALIZE LIFT AND MOMENTS. 

C 

EXLIF-T=EXLI FT/THRUST 
EXPPCH==EXPPCH/(THRUST«0.3 333) 

POVERT=OL!lFT/T.HRUST 
TM=PMOM/ (THRUST*0-3333) 

C 

C CALCULATE WEIGHTED STANDARD DEVIATION. 

C 

TQTFOR=0.0 

totar=o.o 
DO 211 M=l,15 
DO 211 K=5l,25 
IF (M.EQ.15) THET2=l«O.C 

IF(M.NE.lv5) THET2=; (THETA {M+1)+THETA(M) )/2. 

RA2=(«ADIK+ll-«RAD(Kn/2.. 

IF(K.6Q.ir RA1=0.50 


<1 



noon 


IF(M.EQ.l)' THET1=0.0 

IF(K.NE.l): RAI=(RAD(K--1)+RAD<K) )/2. 

If (H.NE. 1); THETl=(THETA(M~l)-»-THETAlM) )/2. 

AREAl=-R*R^(RA2**2-RAi**2 THET2— THETl )*RA01AN/2. 

CPDIFF=P(K,M)-D(M,K) 

IF{A85(CPDIFF).GT.50.I GO TO 211 
FORDIF=(CPDIFF*AREAU=>=«'2 
TOTFOR=TOTFOR+FQRDIF 
T0TAR=ARE=A1=«‘*2+T0TAR 

211 CONTINUE 
STDEV=SQRT(TOTFOR/TOTAR) 

WRITE:(6,212) STDEV 
WRITE(6,3) THRUST 
WRITE{6,2.) TAREA 
WRITE(6fl.y PGVERT,TM 
WRITEI6,A.) EXtlFT^EXPPCH 

212 FORMAT(5X,.*WEIGHTED STANDARD DEVIATION =»,F10.4,//) 

4 FORMAT!///, 5X^» EXPERIMENTAL VALUES* . /5X , •**=»**=i'’f‘*****^=«'*>i'*^^ ’ t g 

1/,5X,*LIFT =*,,F10. 4, /,5X, 'PITCHING MOMENT =',F10,4) 

1 FORMAT I5X',.* TOTAL NON-DIMENSIONALIZED LIFT = • , F10.4 , //, 5X’,^ TOTAL 

1 NCN-DIMEN'SICNALIZED PITCHING MOMENT FIO .4 ) 

2 FQRMAK5X, 'TOTAL AREA =* ,F10. 4, ' SQUARE INCHES',///) 

3 F0RMAT(5X, 'THRUST =»^F10*4) 

RETURN 
END 

SUBROUTINE EXPCP 

SUBROUTINE EXPCP CONVERTS RAW DATA FROM SCANIVALVE PORTS TO VALUES 
CORRESPONDING TO SPECIFIC FIELD POINTS SPECIFIED BY RADIUS AND ANGLE. 

CONMQN/TEN/0121,26) 

DIMENSION CP(10,46) 

DO 3 J=l,21 
DO -3 K=l,26 



non 


3 D<JfK)=100. 

READ EXPERIMENTAL OATA- 

■ READ(5,6) nCPiI»JJ»X^l»10)*j=ltA5),(CPUtA6)*I = lt9) 
6 FaRMATdOFS.At 
OQ 50 J=l,20 
IF(J-GT.16) go to 50 
D( J»1)=CP.U» J) 

D{ Jt2,)=CPClt J+17) 

DlJ,A)=CP5(2,vJ+9 ) 

0U*5)=CPS2t J+30) 

D(J,6)=CP<3, J+1 J 
D< J»7)=CPH3f J+22) 

D( J,9)=CP44» J + 1^ ) 

D{J,ll) = CR<5i.J+6) 

D(J,12)=CPt5»J+23) 

D( J«14)=CR<6,J+15) 

0(J,16)=CP17,J+7) 

0(J,1T)=CP<7,J+28) 

D(J,19)=CPt8,J+20); 

D( J*21)=CRt9tJ+12) 

D( J,22)=CP(9»J+29) 

D( Jf23) = CT>(.10j*J) 

DU»24)=CP( lOrJ + 17) 

IFU.GT.15) GO TO 50 
D( J,10)=CP<4,J+31) 

IF(J.GT.14) GO TO 50 
D^J+7»8)'=CP(4^J) 

DlJ+2*13)^CPt6, J) 

IF(J.6T. 12) GO. TO 50 
D( J+9^20)=:CP^9, J) 

IF(J*GT.10) GO TO 50 
D(J,1S)=CP(6, j+36) 



IF(J..GT.9): GQ TO 50 
D{ J,20)=CP(8f J+37) 
IF<J.GT.8) GO TO 50 
D( J»3)=CPa,J+3fi) 

D( J + 8,3)=CP(>2*J) 
IF(J.GT.7l GO TO 50 
D( JT8).=CR(3,J^'39) 
IFIJ.GT,6); GO TO 50 
0{ J + 15,10) = CP'(5, J) 

0{ J+10,15)=CP(7,J ) 
IFCJ.GT.5) GO TO 50 
D(J + 16,2)’=CP{lf J+33) 
0( J + 16*4)=;CP(2t J+25) 
0{ J+16,6)=CP13, J+17) 

0 1 J+16t 12 )=CP { 5 f J-«;39 ) 
Di J+16fl4l=CP(6,J4.31) 
0{ J+16,16)=CP{7, J+23) 
IF(J.GT.4l GO TO 50 
D(J+3,25)=^CP(10,J+34). 
IF<J«GT.3). GQ TO 50 
Ot J + 8?,25}=,CPilO,J+38l 
50 0( J+l»18)=XPt8t 
D(1,13)=CP(5,45) 
Di2,13)=CP(5,46) 
D(l,18)=CR(7,i46) 
D{21,T)=CP( U17> 
0(21»3)=CP(2,9) 
D(21,5)=CP(3,1) 
D(21,7)=CP(3,39) 
D(21,9)=CP(4,31) 
D(21,11)=CP(5^23) 
D(21*13)=CP16,.15) 
Dt21»15)=CP(7r7 ) 
D(2l»17)=CP(7,.4S) 





0(21, 19)=CP( 8^37) 

D(2I,21)=CP(9^29) 

D(21*22)=CP^9^46) 

D(2l»23)=CP(10,17) 

0(21, 24)-CP( 10,34) 

D(15,25)=CP( 10,42) 

0(5,26)=CP( 10^43) 

D(11,26)=CP( 10,^^) 

D( 15, 26)=CP( 10,45) 

RETURN ^ 

END 

SUBROUTINE WAKE(K,.M) 

C SUBROUTINE WAKE DEFINES TEMPORARY FIELD POINTS FOR USE IN 

C CALCULATING .PRESSURE COEFFICIENTS IN THE WAKE REGION. 

C 

COMMQN/FOUR/XO, YO 
XW=XO 
C 

C DEFINE WAKE BOUNDARY. 

C 

YW= 1 . 819* ( XW.*ffO .9 26 ) 

IF(YO.GT.YW) GO TO 10 
YO=YW 

10 RETURN 
END 
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